Scaling of spontaneous rotation with temperature and plasma current in tokamaks 
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Using tlieoretical arguments, a simple scaling law for the size of the intrinsic rotation observed 
in tokamaks in the absence of momentum injection is found: the velocity generated in the core of a 
tokamak must be proportional to the ion temperature difference in the core divided by the plasma 
current, independent of the size of the device. The constant of proportionality is of the order of 
10km -s"^ • MA • keV~^. When the intrinsic rotation profile is hollow, i.e. it is counter-current 
in the core of the tokamak and co-current in the edge, the scaling law presented in this Letter 
fits the data remarkably well for several tokamaks of vastly different size and heated by different 
mechanisms. 

PACS numbers: 52.25.Fi, 52.30.-q, 52. 55. Fa 



Introduction. Due to their axisymmetry, tokamak plas- 
mas can be made to rotate at high speeds if momen- 
tum is injected into them. If the rotation is sufficiently 
large, large scale magnetohydrodynamic (MHD) instabil- 
ities are stabilized [1,] and the turbulent transport of en- 
ergy can be much reduced 0-01 • Unfortunately, ITER 
[sl, the largest magnetic confinement experiment cur- 
rently being built, is not expected to have effective mo- 
mentum deposition due to its size and high density. As a 
result, there has been mounting interest in the intrinsic, 
or spontaneous, rotation observed in tokamaks without 
momentmn injection Q. If this intrinsic rotation could 
be made large, it could be used to prevent instabilities 
and reduce turbulence as is done with momentum injec- 
tion. Understanding the origin of this rotation is also an 
interesting physics question. This has driven several ex- 
perimental [iHlot , numerical 11 - 1^ and theoretical [l^ 
17[ studies. So far, numerical results for intrinsic rotation 
have only been obtained using global gyrokinetic simu- 
lations that have been recently proven to be flawed for 
radial momentum transport in the core of tokamaks [l8| . 

In this Letter, we use very simple theoretical argu- 
ments to show that the velocity difference within the core 
of a tokamak must scale proportionally to the ion tem- 
perature divided by the plasma current. The constant 
of proportionality is independent of machine size and is 
of order c^/e = 10km • s~^ • MA • keV~^, where c is the 



tSee the Appendix of F. Romanelli et al., Proceedings of the 23rd 
IAEA Fusion Energy Conference 2010, Daejeon, Korea. 



speed of light and e is the proton charge. We show that 
pulses with hollow intrinsic rotation (counter-current ro- 
tation at the magnetic axis and co-current at the edge) 
from machines whose sizes range from tens of centime- 
ters to several meters, that have very different plasma 
currents (from 0.1 MA to 2.5 MA), and that are heated 
by different mechanisms (JET DIII-D fg^, TCABR 
and TCV 7]) follow the theoretical scaling. 

Theoretical arguments. In a tokamak plasma, turbu- 
lence and collisions transport momentum across magnetic 
surfaces. Momentum can be injected with neutral beams 
and radio frequency waves (RF) [i^l, but in many occa- 
sions there is no external source of momentum. When 
the latter is the case, the toroidal angular momentum 
flux n through every flux surface must be zero, even 
though significant rotation can often be observed exper- 
imentally. Only the angular momentum in the toroidal 
direction is relevant. In the poloidal direction, the flow 
is strongly damped by coUisional processes, which pass 
the momentum through the magnets to the structure of 
the tokamak. Thus, to calculate intrinsic rotation pro- 
files, it is necessary to calculate the dependence of 11 on 
the toroidal rotation frequency VL^ and then solve the 
equation 11(17^) = for Vl^. 

Both turbulence and collisions occur on time scales 
that are longer than the inverse of the gyrofrequency, 
which means that the particle trajectories can be under- 
stood as a fast gyromotion around guiding centers, which 
move fast along magnetic field lines and drift slowly 
across them. This is the physical idea underlying gyroki- 
netics, which is the most commonly used approximation 
in transport simulations 21-2^. 
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Even in the absence of turbulence and collisions, par- 
ticles move out of the surface of constant magnetic flux 
where they started due to the VB and curvature drifts, 
but they remain within a given distance of it. This dis- 
tance is of the order of the poloidal gyroradius pe = 
mcvth/eBg, where e and m are the charge and mass of 
the particle, Vth is the thermal speed, and Eg is the 
poloidal component of the magnetic field. Note that 
pe = {B/Bg)p, where p is the particle gyroradius and 
B is the total magnetic field. In most tokamaks, B/Bg is 
of order 10. Tokamaks are constructed so that pg Lt, 
where Lt is the characteristic length of variation of the 
temperature T, which we use as our length of reference. 

Collisions cause transport, known as neoclassical trans- 
port [2^ , because each collision makes the particle move 
from one drift orbit to another separated by pg. Tur- 
bulent transport is caused by electromagnetic fluctua- 
tions, of which the most virulent are believed to be driven 
by the ion temperature gradient (ITG). For ITG turbu- 
lence well above marginality, the characteristic correla- 
tion length is {B / Be){a/ Lt)p ^ pg, where a is the minor 
radius of the tokamak {Lt ~ a). This scaling is not 
based on the drift orbits as is in the case of coUisional 
transport, but on critical balance between the parallel 
and the perpendicular dynamics [2^ . It is observed in 
experiments that the turbulent transport scales approx- 
imately linearly with B/Bg [2^, as predicted by critical 
balance 

In general, tokamaks are geometrically up-down sym- 
metric to a great degree in the core. In such tokamaks, 
to lowest order in pg/Lx, the transport of momentum 
can only be different from zero if a preferred direction is 
given by either rotation or rotational shear. The lowest 
order cancellation of the radial momentum flux in the 
absence of rotation is due to a fundamental symmetry 
of the turbulence and the particle motion 29l431j . Here 
we are assuming il^ vth / R, this being the ordering for 
which the rotation and its shear enter in the lowest or- 
der gyrokinetic equation (3^.[33j. Thus, schematically, to 
lowest order in pg/Lx, 
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where r is the radial coordinate, R is the major radius, 
i^t is the turbulent viscosity, — t't/^pinch is the turbulent 



pinch of momentum [34 , 



and i^c is the collisional vis- 
cosity. Equation ([1]) has the main features of momentum 
transport in up-down symmetric tokamaks: momentum 
transport can only happen when ^ or d^hJdr ^ 0, 
and it changes sign when fi^ and /dr do [3l| . It can 
be thought of as a Taylor expansion of the complicated 
function W{d^l^/ dr, fl^) around fi^ = and dfl^/dr = 0. 

The equation for intrinsic rotation is 11 = 0, and with 
the lowest-order expression ((!]) for 11, the solution is 
il^ oc exp[-/ dr£~.]^^^{l + Uc/vt)~^]- It is then possi- 
ble to obtain intrinsic rotation if rotation is generated in 
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FIG. 1: Sketch of drift orbits. 



some region of the plasma (for example, in the edge) and 
pinched to other regions. However, this mechanism is not 
fully satisfactory because it cannot explain the variety of 
observed proflles In particular, VL^ cannot change 

sign, contradicting experimental observations (as we will 
show in the next section). Unfortunately, to lowest order, 
Eq. ([1} is correct and no other mechanism for intrinsic 
rotation can be obtained. 

If the expansion in pg/Lx ^ 1 is continued to next or- 
der, the rotation and its shear are not the only physical 
factors that provide a preferred direction and can give 
raise to momentum transport: the pressure and temper- 
ature gradients also break the up-down symmetry. Con- 
sider the guiding centers of two particles (1 and 2) that 
at point A at the outboard midplane of a tokamak have 
velocities in opposite directions, as sketched in Fig. [T] 
The dashed line represents the cut of a surface of con- 
stant magnetic flirx through a poloidal plane (the axis of 
symmetry is the dash-dot line). The poloidal magnetic 
field Bg is parallel to the dashed line and points counter- 
clockwise, whereas the toroidal magnetic field -B^ points 
towards the reader. At point A, particle 1 (red orbit) 
travels counterclockwise, and since to lowest order it fol- 
lows the magnetic field, its toroidal velocity v^i is point- 
ing towards the reader. Particle 2 (blue orbit) travels 
in the opposite direction both poloidally and toroidally. 
Orbits separate from the flux surface a small distance 
of order pg. Particle 1 moves towards the center of the 
tokamak because its poloidal velocity is counterclockwise. 
Particle 2 drifts outwards. Because of the temperature 
gradient, the center of the tokamak is hotter, and parti- 
cles like particle 1 will have more energy, of the order of 
{pg / LT)mv'^f^, breaking the symmetry and, in this sim- 
plifled picture, making the plasma rotate counterclock- 
wise poloidally, and towards the reader toroidally. Fig. [T] 
shows that whereas the direction of the magnetic field is 
unimportant, the vector B x VT does give a preferred 
direction at higher order in pg/Lx parallel to or against 
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which the plasma will tend to rotate. The mechanism de- 
scribed here does not determine the sense of the toroidal 
rotation, but it does demonstrate that background gradi- 
ents break the up-down symmetry and that the effects of 
this symmetry breaking are of order pq/Lt- Calculating 
all these effects is a rather sophisticated analytical task, 
involving many factors subtler than the simple argument 
given above 1461 ]. 

The next-order contributions to momentum transport 
in p h/Lt ^ 1 were first calculated in neoclassical theory 
[36l |37| , where they are proportional to radial derivatives 
of the ion temperature. Models to calculate the next- 
order contributions to turbulent transport have also been 
proposed 16|, |17| . Near marginality, the turbulence am- 
plitude is small, and the neoclassical corrections to the 
distribution function of order pg / Lt due to finite drift or- 
bit size are the dominant mechanism that breaks the up- 
down symmetry of the turbulence. Well above marginal- 
ity, the characteristic eddie size is pg [l^l, allowing the 
turbulence to sample regions in which the temperature 
gradient differs by pg/Lx, and breaking the symmetry 
this way. In general, we expect the new next-order terms 
to depend strongly on density and temperature gradi- 
ents because these drive the turbulence. Schematically, 
as shown in we may write 
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From ([3]), setting 11 = and assuming that the 
scale length of ft^ is of order Lt, we obtain Rfl^ ~ 
(pe/LT)vth - {c/eBg){T/LT), where T = mv^j2 is the 
temperature. The poloidal magnetic field is given by the 
toroidal plasma current Ip, Bg ^ Ip/cLs, where Ls is 
the characteristic length of variation of Eg. Therefore, 
Rfl^ - {Lb/Lt){c^ /e){T/Ip). In the core. Lb and Lt 
are both of the order of the minor radius a, so the toroidal 
velocity is 



V4, = RUi 



e /„ 



(3) 



This equation gives the scaling of intrinsic rotation in 
the core with temperature and plasma current. It is in- 
dependent of machine size. The dimensional constant of 
proportionality is c^/e = 10km • s~^ • MA • keV~^. 

Experimental measurements. We now compare exper- 
imental data from different machines that show similar- 
ities in their intrinsic rotation profiles. In Fig. [2] two 
pulses from JET represent two distinct types of intrinsic 
rotation profiles: the ones in which the toroidal velocity 
increases from the magnetic axis towards the edge of the 
tokamak (red profile), which we call hollow profiles, and 
the ones in which it decreases (blue profile), which we call 
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FIG. 2: Intrinsic rotation profiles (a) and ion temperature 
profiles (b) in JET plasmas with ICRH, pulses 66395 (red) 
and 74692 (blue). The rotation in the co-current direction is 
positive rotation. The position of the magnetic axis is around 
ii = 3 m, the separatrix is around R = 3.8 m. 



peaked profiles (the toroidal velocity is deemed positive 
if it is co-current). The two pulses in Fig. [2] have very 
different input power and plasma current, and they are 
only meant to be examples of the two types of velocity 
profiles. The peaked profiles need not have higher tem- 
perature gradients than the hollow profiles. The velocity 
at the edge is mostly co-current, and this seems to be 
common to all tokamaks with low magnetic ripple in the 
absence of momentum injection. In JET, the hollow pro- 
files correspond to Ohmic shots and some of the Ion Cy- 
clotron Resonance Heating (ICRH) pulses in both Low- 
Confinement Mode (L-mode) [10] and High-Confinement 
Mode (H-mode) [s^l- The cases with peaked profiles are 
all ICRH L-mode and H-mode shots. 

To check ([3]), we compare the pulses with hollow core 
velocity profile for four different tokamaks: JET [lo| . 
DIH-D 9], TCABR ^ and TCV 0. To character- 
ize the velocity generated intrinsically in the core, we 
use the difference in toroidal velocity AV^ between the 
minimum of toroidal velocity closest to the magnetic axis 
on the outboard side and the first maximum encountered 
when moving from the magnetic axis towards the edge on 
the outboard side. This definition of AV^ is illustrated 
in Fig. [2ja). The parameter AV^ attempts to exclude 
any intrinsic velocity generated at the edge - most likely 
by means not covered in our theoretical discussion above. 
To give a measure of the sources generating intrinsic rota- 
tion in we use the difference in ion temperature AT^ 
between the magnetic axis and the temperature at the 
top of the pedestal in H-modes, or the temperature mea- 
surement that is the closest to the separatrix in L-modes. 
The difference AT^, illustrated in Fig. [21[b), excludes the 
ion temperature jump in the pedestal in the case of H- 
modes. Fig.[3]shows /pAV^ vs. AT^ for various tokamaks 
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FIG. 3: Toroidal velocity difference in the core AVri) multiplied 
by plasma current Ip against the ion temperature difference 
ATi in the plasma core. The line is the least-square fit of the 
data to @. The slope is 18km ■ s"^ ■ MA ■ keV~V 



[43|. According to ([3|), we expect 



L,AVs = a- AT,,. 



(4) 



The dimensionless prefactor a could not be determined 
in our qualitative theoretical discussion, but we can find 
its value from the present experimental analysis. The 
data is consistent with an approximate linear dependence 
with a slope of (18 ± 4) km • s~^ • MA • keV~^ for all ma- 
chines, giving a ~ 1.8 ± 0.4. The slope was determined 
by least-square fitting and the error is the 99% confidence 
interval. In Fig. [3] there are Ohmic, ICRH and Electron 
Cyclotron Resonance Heating (ECRH) shots, L- and H- 
modes, plasma currents spanning from 0.1 MA to 2.5 
MA, and machines of sizes ranging from tens of centime- 
ters (TCABR) to meters (JET). The fact that both the 
scaling and the prefactor seem to be valid for this variety 
of situations suggests that the theoretical ideas proposed 
above are robust. 

When the same analysis was attempted for the peaked 
profiles in JET, the trend was not so clear, but we cannot 
conclude that the scaling is absent either. The study of 
peaked profiles will be the object of a future publication. 
Here we review several pos sible explanations for the lack 
of a clear scaling. In a change from Trapped Elec- 

tron Mode (TEM) driven turbulence to ITG turbulence 
was proposed as the cause for the transition between 
peaked and hollow profiles. If this is the case, ATi is 
not a good parameter to work with because TEM turbu- 
lence depends strongly on the electron density gradient, 
for which ATi is not a good proxy, and on the electron 
temperature profile that for low coUisionality may differ 
from Ti. It is also possible that the peaked-profile cases 



are dominated by the inward pinch of momentum gen- 
erated at the edge [l^l, making the rotation in the core 
correlated to the parameters at the edge and not to the 
parameters of the core. With our preliminary analysis of 
peaked profiles in JET, we cannot decide if the transi- 
tion from hollow to peaked profiles is due to the reasons 
given above, or other reasons not considered here. Ex- 
periments in DIII-D show that shaping affects intrinsic 
rotation, with high triangularity shots tending to have 
peaked profiles. Shaping may affect ITG and TEM tur- 
bulence differently, making one or the other type of tur- 
bulence dominant for high triangularity and hence decid- 
ing the direction of rotation in this way, or it may have 
an effect on the direction in which ITG or TEM turbu- 
lence pushes the plasma. Further study is needed. Even 
though the trend with ATi and Ip was not so clear, the 
velocity difference AV^ was stiU of the same order as ([3]). 

Discussion. Using simple theoretical arguments, we 
have shown that the intrinsic rotation generated in the 
core must scale according to ([4]). Hollow intrinsic ro- 
tation profiles from very different tokamaks follow this 
scaling. The scatter in Fig. [3] is to be expected since 
(|4|) is derived from an order of magnitude estimate and 
pref actors of order unity may vary from shot to shot. 

There are ways of generating intrinsic rotation that 
have not been considered in this Letter. For exarriple, in 
the core, RF heating can transport momentum 3^ i^l 
due to the large orbits of energetic ions. In the edge, 
direct particle losses can generate rotation [41| . It seems 
that these effects are not important in the cases presented 
in Fig. [3] because these include shots with and without 
energetic ions, and with and without a pedestal. We 
do not know how generic this is. We have introduced 
a dimensionless parameter a = elpAV^/c^ ATi, which 
was of order unity for a variety of regimes and machines 
considered here. It would be very instructive to quantify 
experimentally measured rotation in other cases in terms 
of this parameter. In cases that a is significantly larger 
than unity, the rotation must have external origin, such 
as energetic ions, edge effects or momentum injection. 

The experimental results presented above cannot de- 
termine if the transport of momentum is dominated by 
collisions or turbulence because both have the same scal- 
ing ([3]). Since turbulent viscosity is of the same order as 
the thermal diffusivity 35|, |42|444| , and turbulent trans- 
port usually dominates, we expect the pg/ Lt corrections 
to the turbulent momentum transport to play the domi- 
nant role in driving intrinsic rotation. 
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